a particles and the pasta phase 
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The effects of the a particles in nuclear matter at low densities are investigated within three 
different parametrizations of relativistic models at finite temperature. Both homogeneous and inho- 
mogeneous matter (pasta phase) are described for neutral nuclear matter with fixed proton fractions 
and stellar matter subject to /3-equilibrium and trapped neutrinos. In homogeneous matter, a par- 
ticles are only present at densities below 0.02 fm -3 and their presence decreases with the increase 
of the temperature and, for a fixed temperature, the a particle fraction decreases for smaller proton 
fractions. A repulsive interaction is important to mimic the dissolution of the clusters in homoge- 
neous matter. The effects of the a particles on the pasta structure is very small except close to 
the critical temperatures and / or proton fractions when it may still predict a pasta phase while 
no pasta phase would occur in the absence of light clusters. It is shown that for densities above 
0.01 fm 3 the a particle fraction in the pasta phase is much larger than the a particle fraction in 
homogeneous matter. 
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I. INTRODUCTION 

The pasta phase is a frustrated system F|]-[5| present at 
densities of the order of 0.006 - 0.1 fm -3 [6| in neutral nu- 
clear matter and 0.029 - 0.065 fm -3 @,[1] in ^-equilibrium 
stellar matter, where a competition between the strong 
and the electromagnetic interactions takes place. 

The basic shapes of these complex structures were 
named jl| after well known types of cheese and pasta: 
droplets = meat balls (bubbles = Swiss cheese), rods = 
spaghetti (tubes = penne) and slabs = lasagna, for three, 
two and one dimensions respectively. The pasta phase is 
the ground state configuration if its free energy per par- 
ticle is lower than the corresponding to the homogeneous 
phase at the same density. 

A complete equation of state capable of describing 
matter ranging from very low densities to few times satu- 
ration density and from zero temperature to around 100 
MeV is a necessary step towards the understanding of 
the stellar core collapse, supernova explosion and pro- 
toncutron star evolution. The constitution of the pulsar 
crust plays a definite role in the emission of neutrino and 
gravitational waves. In the inner crust of neutron stars 
(zero temperature, very low proton fraction, matter in 
/3- equilibrium) the pasta phase is expected to be present 
and to coexist with a neutron gas. In supernova (finite 
temperature, proton fraction around 0.3) the pasta phase 
is structured in such a way that there is no neutron gas 
or it is very low in density [9J. 

In previous works @, 0] we have studied the existence 
of the pasta phase at zero and finite temperature within 
different parametrizations of the relativistic non-linear 
Walecka model (NLWM) [ljj, namely NL3 [11], TM1 
[12| and GM3 [13J and the onset of the pasta phase 



with different parametrizations of the density dependent 
hadronic model, namely TW [3], DDH5 [l| and GDFM 
[la . Il7f | respectively. In both works two different meth- 
ods were used: the coexisting phases (CP), both at zero 
and finite temperature and the Thomas-Fermi (TF) ap- 
proximation at zero temperature only. It was found that 
matter in /^-equilibrium presents a small (at zero temper- 
ature or small temperatures) or non-existent pasta struc- 
ture (at finite temperature above ~ 4 MeV) Q. In fact, 
even for Skyrme forces it was shown that depending on 
the model, above 2-4 MeV there is no pasta phase for 
/^-equilibrium matter 18[. If /3-equilibrium is imposed, 



the pasta phase could not be found in a CP calcu- 
lation for the same surface energy parametrization used 
for fixed proton fractions. This indicates the necessity to 
use a good parametrization for the surface energy which 
is temperature, proton fraction and geometry dependent, 
as also stressed in jig, [l7| • The specific problem of an ap- 
propriate parametrization for the surface energy is tack- 
led in the present work. 

In the above calculations we have assumed predefined 
shapes for the pasta clusters and considered that the 
groundstate configuration is the one that minimizes the 
free energy. In [19] the authors could show by using ab- 
initio numerical simulations that in fact pasta phases can 
be formed in collapsing supernova. By compressing a bec 
lattice of spherical nuclei it was shown that an ordered 
structure of rod-like nuclei could be formed. 

The importance of the a particles has been pointed 
out in the recent literature |20l - [24| . It is the most 
strongly bound system among all light systems and it 
certainly plays a role in nuclear matter. Lattimer and 
Swesty worked out the equation of state (EOS) in the 
compressible extended liquid drop model based on a 
non-relativistic framework [2(| appropriate for supernova 
simulations, for a wide range of densities, proton frac- 
tions and temperatures, including the contribution of 
ev-particle clusters. An excluded volume prescription is 
used to model the dissolution of a particles at high den- 



sities. The same is done by Shen et al. in [2l|, [25[ 
where non-uniform matter composed of protons, neu- 
trons, a particles and single species of heavy nuclei is de- 
scribed with the Thomas-Fermi approximation and the 
TM1 paramctrization of the non-linear Walecka model 
(NLWM). At low densities, these particles are approxi- 
mated by classical gases. In [26] nuclear statistical equi- 
librium equations are calculated and the a particles are 
also taken into account. In [22| an effective four-body 
equation that includes self-energy corrections and Pauli 
blocking is done in a consistent way. The influence of 
cluster formation on nuclear matter EOS is calculated 
and the occurrence of instabilities investigated in [271 ]. 
In J23| the virial equation of state of low-density nuclear 
matter composed of neutrons, protons and a-particles 
is presented and it is shown that the predicted a parti- 
cle concentrations differ from the predictions of the EOS 
proposed in [2(| and (25(. The inclusion of small clus- 
ters in the EOS is revisited in [2J], where the most im- 
portant thermodynamical quantities are calculated with 
light clusters up to the a particle with a density depen- 
dent relativistic model. The conditions for the liquid- 
gas phase transition are obtained and it is seen how the 
binodal section is affected by the inclusion of these clus- 
ters. Moreover, an EOS is obtained from the gaseous low 
densities with a clusterised matter to higher cluster-free 
homogeneous matter. 

In the present paper we investigate the influence of 
the a particles both on the homogeneous matter and 
on the onset, size of the clusters and structure of the 
pasta phase within the relativistic mean field approxi- 
mation. As all results so far show model dependence, 
three different parametrizations are employed. We have 
chosen the NL3 
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GM1 [13j parametrizations of the 
NLWM and the TW parametrization of the density de- 
pendent hadronic model [lj]. The calculations are car- 
ried out for fixed proton fractions and also for matter in 
/3-equilibrium with trapped neutrinos. Although in pre- 
vious works it was found that matter in /3-equilibrium 
presents a small (at zero temperature) or non-existent 
pasta structure (at finite temperature) that is very sensi- 
tive to the surface energy [y, Q, once trapped neutri- 
nos are present, the picture changes considerably due 
to the large fraction of protons. According to studies 
on binodals and spinodals underlying the conditions for 
phase coexistence and phase transitions [l8l . I28J . non- 
homogeneous matter with trapped neutrinos is expected 
to be found until temperatures around T = 12 MeV, de- 
pending on the model considered. For this reason, we 
investigate this possibility next. 

The paper is organized as follows: in section II we 
briefly review the formalism underlying the homogeneous 
neutral npe matter with the inclusion of the a particles. 
In section III the pasta phase is built with the help of 
the coexisting phases method and the prescription for 
the introduction of the as is given. A complete study for 
the parametrization of the surface energy based on the 
Thomas- Fermi calculation is performed and presented in 



section IV. In section V our results are displayed and 
commented and in section VI our conclusions are drawn. 



II. THE FORMALISM 

We consider a system of protons and neutrons with 
mass M interacting with and through an isoscalar-scalar 
field <fi with mass m s , a isoscalar- vector field V^ with 
mass m v , an isovector- vector field b M with mass m p . We 
also include a particles as bosons with mass M a and 
a system of electrons with mass m e . We do not con- 
sider models with 6 mesons next, but their introduction 
is straightforward. The Lagrangian density reads: 
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the a particles are described as in [24[ by 



£a = r(i^ a )*(»2)^) " ,C««i. 



with 



where 



iD^ 
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(2) 

(3) 
(4) 

(5) 

(6) 

(7) 

and r Q is included for mimicking the excluded volume 
effect and consequent a particles dissolution at high den- 
sities. In [24| a more complete effective interaction for 
the a particles is considered, including the a meson-a 
particle interaction and the Pauli shifts for the binding 
energy. In our approach we calculate a lower bound for 
the a particles present in the medium. As a reference we 
also consider the opposite limit and take the a-particles 
as a gas of free particles. In most of our calculations T a 
is either set to zero or to T a — 4T V . We have also made 
tests using a lower value for this interaction (2r„) and 
the results arc commented when appropriate. 
The electron Lagrangian density is given by 

and the meson Lagrangian densities are 



M„ = AM - B n 
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where £1 M „ = <9 M V„ - d v V^ , and B M „ = <9 M b„ - d^b^ - 
r p (b M x b„). The parameters of the models are: the 
nucleoli mass M = 939 MeV, the coupling parameters 
r s , r„, r p of the mesons to the nucleons, the elec- 
tron mass m e and the electromagnetic coupling constant 
e = y/4ir/137. In the above Lagrangian density r is 
the isospin operator. When density dependent models 
are used, the non-linear terms are not present and hence 
k = A = £ = and the density dependent parameters 
are chosen as in [lj, [jj| [29| ■ When models with constant 
couplings are used, I\ is replaced by gi, where i — s,v,p 
as in [111 . 1 131 ] . The bulk nuclear matter properties of 
the models we use in the present paper are given in Ta- 
ble HJ We also include in the table some properties at 
the thermodynamical spinodal surface: p sp in is the up- 
per border density at the spinodal surface for symmetric 
matter (it defines the density for which the incompress- 
ibility is zero), pt and Pt are, respectively, the density 
and the pressure at the crossing between the cold /3- 
cquilibrium equation of state and the spinodal surface. 
They give a rough estimate of the density and pressure 
at the crust-core transition 6, 8] . For symmetric matter 
or large proton fractions we expect that GM1 will have 
the largest extension of the pasta phase because it has 
the largest value of p sp i n - For cold /3-equilibrium matter 
TW predicts the largest pasta extension, with the largest 
transition density pt- 



TABLE I: Nuclear matter properties at the saturation density 
and at the spinodal surface 





NL3 


GM1 


TW 




[11] 


[13] 


[14] 


B/A (MeV) 


16.3 


16.3 


16.3 


pa (fin -3 ) 


0.148 


0.153 


0.153 


K (MeV) 


272 


300 


240 


Esym. (MeV) 


37.4 


32.5 


32.8 


M*/M 


0.60 


0.70 


0.55 


L (MeV) 


118.3 


93.8 


55.3 


K sym (MeV) 


100.5 


17.9 


-124.7 


Qo (MeV) 


203 


-216 


-540 


K T (MeV) 


-698 


-478 


-332 


p spin (fm~ 3 ) 


0.096 


0.100 


0.096 


Pt (fm" 3 ) 


0.0646 


0.0743 


0.0850 


Pt (MeV/fm 3 ) 


0.396 


0.382 


0.455 



From de Euler-Lagrange formalism we obtain coupled 
differential equations for the scalar, vector, isovector- 
scalar and nucleon fields. In the static case there are no 
currents and the spatial vector components are zero. In 
@ a complete description of the mean-field and Thomas- 
Fermi approximations applied to different parametriza- 
tions of the NLWM are given and we do not repeat them 
here. The equations of motion for the fields are obtained 
and solved self-consistently and they can be read off 6, 8] . 



The above mentioned equations of motion depend on the 
the equilibrium densities p = p p + p n , p% = p p — p n and 
Ps = Ps P + Psn where the proton/neutron densities are 
given by 



Pi 



p 2 dp(f t+ -ft-), i = p,n 



and the corresponding scalar density by 
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l + exp[(e*(p)T^)/T] 



where e* = \/p 2 + M* 2 , 

M*=M-T S 4, , 
and the effective chemical potentials are 

r 
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(9) 



(10) 



(11) 



(12) 



(13) 



T3i = ±1 is the isospin projection for the protons and 
neutrons respectively. The density dependent models in 
the mean field approximation contain a rearrangement 
term T,Q t [15| given by: 

In the description of the equations of state of a system, 
the required quantities are the baryonic density, energy 
density, pressure and free energy, explicitly written in 
[0, Q. We refer next only to some of the quantities. 
The free energy density is defined as: 

T = £-TS, (14) 

with the entropy density : 

S = —{£ + P - PpPp- PnPn)- (15) 

In the present work the a particles are included as 
bosons and their chemical potential is obtained from the 
proton and neutron chemical potentials as in |24l . |26| ] : 



p a = 2 (p p + p n ). 



(16) 



The inclusion of the a particles also give rise to a re- 
arrangement term given by 






^ Pa V , 
dp 



(17) 



which is present in the effective chemical potential of the 
a particles: 



/'<- 



T a V 



(18) 



The density of a particles is 



1 

p* = — 



P 2 dp{f a + -fa-), 



with the boson distribution function given by 

f a ± = - 

-1 + cxp[(e Q =Fi/ a )/T]' 



(19) 



(20) 



where e a = ^/p 2 + M 2 . The a free energy density reads 

•Fa = £ a — TS a , (21) 

where £ a and S a stand respectively for bosonic energy 
and entropy density. 

It is worth emphasizing that a term proportional to 
TiQ a p a is present in the pressure. When non linear 
Walecka type models are used there is no rearrangement 
term and the a effective chemical potential is also given 
by eq. (|16p . For the temperatures we consider in the 
present work no alpha particle condensation occurs, and, 
therefore, the above equations do not contain the con- 
densate contribution. 

As for the electrons, their density and distribution 
functions read: 



Pe = ~2~ / P 2d P(fe+ ~ fe-), 



with 



/e±(r,p,i) = 



1 + exp[(e e =F fJ-e)/T] ' 



(22) 



(23) 



where p e is the electron chemical potential and e e = 
y/p 2 + m 2 . We always consider neutral matter and there- 
fore the electron density is equal to the charge density of 
the positive charged particles (protons and alphas) . 
The free energy density of the electrons reads 



Te 



TS e , 



with 



<$e = ~(£e + Pe- PePe)- 



(24) 



(25) 



III. COEXISTING PHASES 

Two possibilities are discussed next: nuclear matter 
with fixed proton fraction and stellar matter with (5- 
equilibrium and trapped neutrinos. 



A. Nuclear pasta 

As in @, Q, for a given total density p and proton 
fraction, now defined with the inclusion of the protons 
present inside the a particles, the pasta structures are 
built with different geometrical forms in a background 



nucleon gas. This is achieved by calculating from the 
Gibbs' conditions the density and the proton fraction of 
the pasta and of the background gas, so that in the whole 
we have to solve simultaneously the following six equa- 
tions: 

P*(yLy n , M* 1 ) = P n { v " v», M* n ), (26) 
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(27) 
(28) 

(29) 



f(p' p + 2pi) + (1 - /)(^ + 2p I J) = Y pP , (30) 

where I and II label each of the phases, / is the volume 
fraction of phase I: 



,_ P-P" 



(31) 



P - P 
where the total baryonic density is 

p = p p + p n +4p a , (32) 

and Y p is the global proton fraction given by 

P P + 2/? 



Y p = 



(33) 



The density of electrons is uniform and taken as p e — 
Y p p. The total pressure is given by P = P 1 + P e + P a . 
The total energy density of the system is given by 

E = f£' + (1 - f)E" +£e+ Ssurf + Seoul, (34) 

where, by minimizing the sum E sur t +£coui with respect 
to the size of the droplet /bubble, rod/tube or slab we get 
5] £ sur f = 2£ C oui, and 



2F 



£coui = ^W(eW /3 (aD( Pp - p»)) 
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(35) 



where F — f for droplets and F = 1 — / for bubbles, a is 
the surface energy coefficient, D is the dimension of the 
system. For droplets, rods and slabs, 
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(36) 



and for bubbles and tubes the above expressions are valid 
with / replaced by 1 — /. 

Each structure is considered to be in the center of 
a charge neutral Wigner-Seitz cell constituted by neu- 
trons, protons and leptons [21]. The Wigner-Seitz cell is 
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FIG. 1: (Color online) Surface energy as a function of the 
proton fraction for GM1, NL3 and TW- 



a sphcrc/cilindcr/slab whose volume is the same as the 
unit BCC cell. In [2l| the internal structures are associ- 
ated with heavy nuclei. Hence, the radius of the droplet 
(rod, slab) and of the Wigner-Seitz cell are respectively 
given by: 



Rn = 



aD 



47re 2 (p£ + 2pi - ft - 2p'J)^ 



1/3 



Rw = 



Rd 



{i-jy/D- 



(37) 



(38) 



In Fig. [T]the surface energy is plotted as a function of 
the proton fraction for the three models under study. It 
is seen that the models differ a lot. 



B. Surface energy 

The authors of [301 ] have studied how the uncertain- 
ties on the surface energy affect the appearance of non- 
spherical pasta structures. In particular, they have 
shown that for typical values of the surface energy non- 
spherical clusters may occur below the transition density 
to uniform matter. Also the authors of Q state that 
the appearance of the pasta phase essentially depends 
on the value of the surface tension. We have fixed the 
surface tension at different values and confirmed their 
claim in [a, Q ; where the surface energy coefficient was 
parametrized in terms of the proton fraction according 
to the functional proposed in [311 ] . obtained by fitting 
Thomas-Fermi and Hartree-Fock numerical values with 
a Skyrme force. 

In the present paper we have considered this point 
in a more systematic and consistent way. Hence, we 
have parametrized the surface energy according to our 
Thomas- Fermi calculations for the relativistic models un- 
der investigation. We have used the Gibbs prescription 
to obtain the a coefficient which, as extensively discussed 



in the literature [32|, |33j , is the appropriate surface ten- 
sion coefficient to be used. We have obtained for the 
three parametrizations used in this work, namely NL3, 
TW and GM1 parameter sets, the corresponding fittings 
for a. Here, we briefly discuss the method that we have 
employed for fitting a and more specific details about 
the numerical algorithm can be found in [6j. First, we 
obtain a density profile for the two component system 
consisting of protons and neutrons in a semi-infinite one- 
dimensional system, where the axis perpendicular to the 
interface has been defined as the z-axis. In our approach, 
we assume the system inside a large box of radius R so 
that the densities at the left and right boundaries corre- 
spond to asymptotic densities that approach the values 
of uniform nucleon matter (nucleus) in equilibrium with 
a gas of drip particles and a uniform gas of drip particles 
respectively. Hence, we need a recipe in order to extract 
the surface energy from the bulk one. The surface tension 
may be expressed as follows: 



dz[e(z) - e d - e ref (z)} 



(39) 



as discussed in [32|. The quantity e re f(z) corresponds to 
a reference energy density associated to the bulk contri- 
bution, e d — lim z _ i . 00 e(z) and e(z) is the energy density. 
If one considers the Gibbs phase coexistence conditions 
and usual thermodynamic relations one obtains for the 
(Gibbs) surface tension coefficient the expression: 



dz (e(z) - e d - fi n [p n {z) - p nd ] 



~Mp [Pp{z) ~ Ppd]) , 

where p nd = lim2_ > , 00 p n (z) and an analogous definition 
for p pd . An alternative and completely equivalent expres- 
sion for the calculation of a, known as thin-wall approx- 
imation 34-36] is given by the following expression: 




dr 




(40) 



In our parametrization of a the above expression has been 
used since it is more convenient for our Thomas-Fermi 
numerical algorithm 6]. In particular, it is specially ad- 
equate for the parametrization of the dependence of a 
on the temperature. We do not try to parametrize the a 
through standard dimensionless functions which, as dis- 
cussed in |33l ] . do not give good fits for relativistic mean 
field models, specially, for the region of small a values 
which are important for the pasta study. So, we follow 
a more pragmatic way adopting a mathematical formula 
for a that gives accurate results for a broad range of neu- 
tron excess and for temperatures up to lOMeV, which are 
adequate for the studies addressed in the present work. 
The following functional for the surface tension coeffi- 
cient, a, is used, 



a = a(x, T = 0) [1 - a(T) xT - b(T)T 2 ] 



(41) 



where x — S 2 stands for the squared relative neutron 
excess: 



5 = P^—PlL =1-2Y P , 



where p, p n , and p p are defined at z = — oo. A table 
with the coefficients a(x,T — 0), a(T) and b(T) for the 
NL3, TW and GM1 parametrizations is shown in the 
appendix. The proton fraction considered throughout 
the calculation of a is the one of the denser phase. 



C. Stellar pasta for matter in /3 equilibrium with 
trapped neutrinos 

In this case, hadronic matter is in /3 equilibrium and 
the electron neutrinos are trapped. The condition of (3 
equilibrium in a system of protons, neutrons, electrons 
and neutrinos is 



M P = Pn - Me + AV • 

Besides the Gibbs conditions given in eqs. 
and P5|) we also impose 



Y 



Pv c 



P 



(42) 

, urn 

(43) 



where Yj is the lepton fraction and has been fixed as 0.4 
and p Vc is the density of electron neutrinos given by 



Pu c = 2~2 / P 2d P(f"e+ ~ Ue-)> 



with 



f«.± = 



l + cxp[(pTAV)AT 
The neutrino free energy density reads 



(44) 



(45) 



(46) 



where £„ c and S Vs stand respectively for the neutrino 
energy and entropy density. 

The densities of interest to the study of the pasta phase 
are too low for the muons to appear, which generally 
occur for densities above 0.1 fm -3 [18| and hence are not 
considered in the present work. 



IV. RESULTS AND DISCUSSIONS 

In this section we show results for the fraction of a- 
particles in homogeneous matter and pasta-like matter 
and discuss the effect of the a-particles on the structure 
of the pasta phase. 
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FIG. 2: (Color online) a fractions as a function of density for 
a) symmetric matter and several temperatures and b) T=5 
MeV and several proton fractions, obtained with NL3. The 
thin curves are for free a-particles and the thick ones include 
the aj-meson-a particle interaction. 



A. a particles in homogeneous matter 

We show the amount of a particles present in homoge- 
neous nuclear matter described by the NL3 parametriza- 
tion for symmetric matter and different temperatures in 
FigEK and for T — 5 MeV and three proton fractions, 
Y p = 0.5, 0.3, 0.1 in Fig. [5b- In all figures the curves 
drawn with thin lines were obtained for free a-particles, 
T a = and the curves with thick lines for T a — 4T„, i.e., 
with the inclusion of the term that mimics the excluded 
volume factor. As mentioned before, this term dissolves 
the clusters. In both calculations the a-particle vacuum 
mass was used. The maximum of the distribution occurs 
for the density that maximizes 2(p p +/x n ) — T a Vo — M a . 
The curves obtained with this term show the same be- 
havior found in Fig. 15 of [2J|, i.e., the a particle frac- 
tion decreases with the increase of the temperature. One 
can see from FigJ^b that the total proton fraction Y p 
has a different effect on the a particle distribution. For 
T a — 4F V , the maximum of the distribution occurs at the 
same density for all fractions, although the amount of a 
particles shows a lower value for the smaller value of Y p . 
Moreover, the density of dissolution of a particles does 



not seem to be sensitive to the proton fraction, because 
it is denned by the isoscalar vector interaction. The a 
particle frations for free a particles in Fig. [2^. and 
[2Jd show discontinuities for large proton fractions 
at T=5 MeV, just above p = 0.02 fm 3 . For lower 
proton fractions no discontinuity occurs. The re- 
gion where the discontinuity occurs is a region 
of instability since the free energy has a negative 
concavity. For T=5 MeV and y p =0.3 there is a 
smooth transition from a region with a large al- 
pha particle fraction to a region with low alpha 
particle fraction. If we increase the proton frac- 
tion the transition is not any more smooth and we 
obtain a jump in the a particle fraction. The dis- 
continuity occurs when the alpha particle fraction 
becomes close to one and the chemical potential 
of the a-particle close to the alpha particle mass. 
For larger densities the nucleon chemical poten- 
tial decreases due to the attraction induced by the 
a meson. The minimum of the nucleon chemical 
potential at p ~ 0.1 fm 3 corresponds to the min- 
imum of the alpha particle fraction. The discon- 
tinuity or kink occurs at the density, or slightly 
above, for which the nucleon chemical potential 
in the presence of alpha particles becomes equal 
to the nucleon chemical potential without alpha 
particles. To correctly treat this region of den- 
sities, where the free energy has a negative con- 
cavity, we should consider that matter is not ho- 
mogeneous and there are two phases. However, 
describing alpha particles as free particles is not 
realistic and when the interaction between parti- 
cles is included the discontinuity in the a particle 
fraction does not occur. 

In Fig. [3^ the amount of a particles present in mat- 
ter with /3-equilibrium and trapped neutrinos with the 
NL3 parametrization is displayed for different tempera- 
tures. Just like in Fig. [2^,), it is seen that the larger 
the temperature the lower the maximum of the a parti- 
cle distribution. With the increase of temperature, the 
maximum is shifted to larger densities as well as the dis- 
solution of the clusters. This behavior was also obtained 
in [241 ]. However in [24[ other smaller clusters, besides the 
a particles, have also been considered. In Fig. |3Jd) the 
same is shown for T=5 MeV and the three models used 
in the present work. TW is responsible for the smallest 
amount of a particles, whilst GM1 for the largest. These 
behaviors are defined by the model properties at subsat- 
uration densities. This can be seen in Fig. [3};) where the 
nucleon effective mass of NL3, GM1 and TW are plotted 
for subsaturation densities. If the nucleon mass is larger 
it becomes energetically more favorable to produce more 
alpha particles and less free nucleons. Both the cluster 
dissolution and the particle maximum occur at smaller 
densities for TW and larger for GM1 for the same reason. 

The virial expansion gives a model independent predic- 
tion of the low density limit of the EOS [23]. In Fig. @]we 
compare the a particle fraction obtained with the virial 
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FIG. 3: (Color online) a fractions as a function of density 
obtained with a) model NL3 for different temperatures; b) 
models NL3, GM1 and TW for T = 5 MeV for matter in 
/3-equilibrium with trapped neutrinos. The thin curves are 
for free a-particles and the thick ones include the u;-meson-« 
particle interaction. In c) the nucleon effective mass within 
NL3, GM1 and TW is plotted at subsaturation densities. 



EOS of low density nuclear matter [23j with the predic- 
tion of the models GM1, NL3 and TW with interacting a 
particles, T a — AT V . We also include the NL3 results for 
free a particles (thin red line). The GM1 parametriza- 
tion with interacting a particles gives the closest results 
to the virial expansion and TW deviates the most. 
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FIG. 4: (Color online) a fractions as a function of density 
obtained with GM1, NL3 and TW models are compared with 
the virial EOS results of ref. [H. 
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FIG. 5: (Color online) a fractions as a function of density 
obtained with model NL3 for a) Y v - 0.3 and T = 5 MeV; 
b) for Y p — 0.3 and T — 8 MeV. The thin curves are for free 
a-particles and the thick ones include the a;-meson-a particle 
interaction. 



B. a particles in the pasta phase 

In the sequel, we analyse the effect of the pasta phase 
on the fraction of a-particles. The a particles are present 
in the background gas which surrounds the pasta struc- 
tures. 



FIG. 6: (Color online) a fractions as a function of density 
obtained with the models under investigation for Y p — 0.5 
and Y p = 0.3 and T = 5, 8 MeV. The thin curves are for T=5 
MeV and the thick curves for T=8 MeV 



Q. 



a 

Q. 



10 

10" 

10" 

10" 

10" 

10" 

10" 

10" 
10" 

10" 
10" 



10 



12 - 



4 - 



10 



12 - 



14 



" ' _^^^ ~^N\ 


" 


NLS^. 


pasta 


T= 5 MeV V\ 


- 


Y 


\ a) 


■Yp-0.5- \ 




Y p = 0.3 — hom.V, 




/^NL3 ^\\ 


pasta 


iy^ T= 8 MeV \ N 






\ b) 


hom.\ 


\ 



10 10 

P(fm" 3 ) 



10" 



FIG. 7: (Color online) a particle density for a) T=5 MeV and 
b) T= 8 MeV and Y p = 0.3 and 0.5 obtained with NL3 for 
homogeneous matter (thin lines) and gas phase of the pasta- 
like matter (thick lines). All calculations include the o;-meson- 
a particle interaction. 



In Fig [5] we have plotted the a particle fraction Y a 
for T=5 and 8 MeV with the proton fraction Y p = 0.3. 
Notice that Y a — 4p Q /p is the fraction of nucleons clus- 
tered in a-particles. In both figures we include the homo- 
geneous matter results (dashed lines) and the pasta-like 
matter results (full lines). We have considered both free 
a particles (thin lines) and interacting a particles (thick 
lines) . If free a particles are considered the fraction Y a in 
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FIG. 8: (Color online) a particle density in the gas phase 
of pasta-like matter for T=5 MeV and /3-equilibrium matter 
with trapped neutrinos obtained with NL3, GM1 and TW. 
The a particle density is given for three cases: free a particles 
T a — (thin lines), interacting a-particles with T Q = 4X„ 
(thick lines), and just for NL3, interacting a-particles with 
T a — 2r„ (red medium thick line). 



homogeneous matter becomes quite high and it may be 
close to one which is quite unrealistic. In pasta-like mat- 
ter the maximum of Y a is smaller, but still quite high: 
it varies between ~ 0.01 and 0.1 and, for Y p = 0.3, it 
is larger for larger temperatures, T = 8 MeV. This last 
tendency is still true when we include the interaction of a 
particles with the w meson. If we analyze the results in- 
cluding the a — uj interaction some important conclusions 
can be drawn: while for homogeneous matter Y a < 10~ 12 
for p > 0.02 fm~ 3 , for pasta-like matter this only occurs 
at the crust-core boundary for p ~ 0.1 fm~ 3 . 

We may consider the prediction of Fig. [5]a lower bound 
prediction since for the a particle interaction we have 
only considered a repulsive interaction which mimics the 
volume exclusion effect. If attraction would have also 
been considered we would expect larger fractions [24J . 
Also we have taken for the coupling constant the nucleon- 
meson coupling constant multiplied by the mass number 
and this could be a too strong coupling. In pasta-like 
matter for symmetric matter the a particle-fraction does 
not change due to the method used to calculate the pasta 
phase. Phase equilibrium for symmetric nuclear mat- 
ter predicts the same proton fraction, equal to 0.5, as 
well as the same pressure for both the dense (liquid) and 
gas phases of the pasta. We may expect that within a 
Thomas- Fermi calculation the a-particle fraction changes 
with density. 

From Fig. [5] it is seen again that the a-particle frac- 
tion is model dependent: GM1 predicts generally larger 
fractions although the difference depends on the density. 
The differences can be as high as one or two orders of 
magnitude. The magnitude of Y a depends on the frac- 
tion of protons present in the background gas. The low 
density onset of the pasta-like matter in Figs. 5 and 6 
lies slightly below the homogeneous matter curves. The 
difference is smaller for the larger temperatures due to 



the distillation effect. The pasta-phase starts at low den- 
sity with small droplets of dense matter in a background 
gas. As soon as the droplets set on the distillation effect 
makes matter in the droplet more proton rich and in the 
gas more proton poor. Therefore, fixing the global pro- 
ton fraction, we expect that at the onset of the droplets 
(the volume fraction of the dense part is very small) the 
gas phase has a smaller proton fraction than the homoge- 
neous matter at the corresponding density. This effect is 
smaller for the larger temperatures and can be observed 
from Figs. 5 and 6. 

The a particle densities are plotted in Fig. [7] for two 
proton fractions (Y p = 0.5 and 0.3) and two temperatures 
(T = 5, 8 MeV). We include the calculation for both 
homogeneous matter and pasta-like matter and we only 
consider the case of interacting a particles. This figure 
gives a hint on the possible effects of a particles in the 
inner crust of a compact star. The effect is larger for the 
larger temperatures and larger proton fractions. Due to 
the existence of a non-homogeneous phase the presence 
of the a-particles extends to larger densities. In fact 
the a clusters dissolve at p ~ 0.01 fm -3 in homogeneous 
matter. If we consider the pasta phase the background 
gas does not exceed this value of the density and, there- 
fore, there are still alpha-clusters at much larger densities 
because they only exist in the gas phase. However, we 
should point out, that with the strongly repulsive inter- 
action considered, Y a takes very small values. 

In Fig [2^.) and [3^) it is clearly seen that the 
maximum fraction of alpha particles decreases 
with temperature, however clusters dissolve at 
larger densities and above 0.01 fm~ 3 the fraction 
of a particles is much larger for T=10 MeV than 
for T=5 MeV. This effect is still seen in the pasta 
phase: above 0.01 fm 3 a particle fractions are 
larger at T = 8 MeV than at T = 5 MeV. 

It is worth emphasing some points related to the effects 
of the a particles present in the inner crust of a com- 
pact star. They are small or negligible in calculations 
of the evolution of protoneutron stars due to the very 
low densities involved. However, simulations of super- 
novae explosions seem to indicate the the internal struc- 
ture of the pasta phase plays a decisive role in avoiding 
the stalling process. Moreover, at the inner crust the 
shear viscosity, thermal conductivity and electrical con- 
ductivity are mainly determined by electron-ion scatter- 
ing. Above neutron drip it is important also to consider 
electron- neutron scattering [37H39I |. It has been shown 
[40J that the shear viscosity in the inner crust, which 
is mainly determined by electrons, is affected by the 
electron-impurity scattering at low temperatures. Re- 
cently, taking into account electron-proton scattering it 
was shown in [39f that the pasta structures would not 
increase as expected the shear viscosity and thermal con- 
ductivity due to an effective reduction of the proton num- 
ber of the clusters owing to ion screening effects. We ex- 
pect that the presence of alpha-particles will give an ex- 
tra contribution that should be calculated. In this paper 
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TABLE II: Densities of the inner edge of the crust (crossing points) for NL3 







EOS 
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fin" 3 


pi J p2 

MeV/fm 3 


pasta structure 
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= 5 MeV 












Y P 
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(*) 


0.004 / 0.077 0.039 / 2.065 


droplets, rods, slabs 
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J- va — 





0.010 / 0.077 0.155 / 2.069 


droplets, rods, slabs 






r - 
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0.004 / 0.077 0.036 / 2.065 
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Y P 


= 0.3 
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(*) 
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0.006 / 0.082 


0.044 / 1.138 


droplets, rods, slabs 






r 

J- va — 


4r„ 
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it was also shown that the scattering of neutrinos from 
neutrons would be defined by the difference between the 
neutron content of the pasta structures and the dripped 
neutron gas. The reduction of the background neutron 
gas due to the alpha particle formation could increase the 
effective neutron content of the pasta structures. 

In order to show the effect of the a particle interaction 
on the background gas of the pasta phase, we plot in Fig. 
IH1 for matter with trapped neutrinos and a lepton fraction 
Yi = 0.4, the density p a obtained with free and inter- 
acting particles. The a-particle density is lower for TW. 
A smaller a-particle fraction is due to the smaller sym- 
metry energy for TW at subsaturation densities, which 
does not favor the formation of isospin symmetric clus- 
ters. GM1, on the other hand, has the largest symmetry 
energy and predicts much larger a-particle fractions. We 
also conclude that knowing with accuracy the strength 
of the interaction is crucial to obtain the correct order 
of magnitude for Y a . At p = 0.08 fin -3 , which corre- 
sponds approximately to the crust-core transition, the 
difference between a free a particle gas and an interact- 
ing one with T Q = 4T V could be as high as 6 orders of 
magnitude. Reducing the repulsive w — a interaction to 
one half, T a = 2T V , reduces the difference to 3 orders of 
magnitude. This test is shown only for NL3 and is rep- 
resented in the figure by the red medium thick line, close 



to the GM1 result with T a = 4T„. 



C. Effect of the a particles on the phase transitions 

In Fig. |H]we have plotted the free energy for homoge- 
neous and pasta like matter described by the NL3 model 
with T=5 MeV and proton fractions Y p — 0.5,0.3 and 
for /3-equilibrium matter with trapped neutrinos for a 
lepton fraction Yi = 0.4. The extension of the pasta free 
energy is determined by the binodal surface and becomes 
smaller when the proton fraction decreases. The system 
in equilibrium chooses the configuration with lower free 
energy so the pasta-like matter defines the ground state 
of the system only while its free energy is lower than the 
free energy of the corresponding homogeneous matter. 

In Tables [Til an d IIIII the densities at the crossing be- 
tween the homogeneous and the pasta phases are given 
for NL3 and TW. The two crossing phases correspond to 
homogeneous matter to pasta {p\ in the Tables), at very 
low densities, and pasta to homogeneous matter (p 2 in 
the Tables), at slightly higher densities. The correspond- 
ing pressures (P 1 and P 2 obtained in the pasta phase) 
are also given. The (*) means that the numbers shown 
for the crossing densities with NL3 and no a particles are 
different from the ones shown in [6| because a different 
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FIG. 9: (Color online) Free energy for the homogeneous and 
pasta like matter obtained within the NL3 for T=5 MeV and 
several proton fractions. The red dashed line defines the pasta 
free energy and its extension is determined by the binodal. 
The crossing of the pasta and the homogeneous free energy 
define the extension of the pasta phase. 



parametrization for the surface energy was used here. 

It is seen from Tables [XT] and Mil that the transition 
density for symmetric matter is smaller than the transi- 
tion density for matter with Y p = 0.3 while the opposite 
would be expected since the binodal has a smaller exten- 
sion for the smaller proton fractions. In fact a Thomas 
Fermi calculation predicts respectively 0.104 and 0.092 
fm -3 for Yp = 0.5 and Y p = 0.3 for NL3 and T= 5 MeV 



41| . As can be seen from Fig. [9] there is a range of den- 
sities close to the transition density when the free energy 
of pasta and homogeneous matter do not differ much. 
The same occurs with a Thomas Fermi calculation. This 
corresponds to the range of densities which is sensitive to 
the method used and the self-consistency of the method. 
One large drawback of the coexisting phases method is 
the fact the proton distribution inside the pasta structure 



is not perturbed by the Coulomb field. In particular, for 
symmetric matter the proton fraction of the structure is 
0.5, larger than the prediction of a Thomas Fermi cal- 
culation which would predict that protons spread giving 
rise to a smaller proton fraction, since the neutrons are 
not affected. This effect is larger the larger the proton 
fraction and is particularly critical for symmetric matter. 
We have verified that the density transition for symmet- 
ric matter at T=5 MeV would increase to 0.087 fm~ 3 if 
the central proton fraction of the structure would reduce 
by 20% due to the Coulomb force. 

The presence of the a particles has only a small ef- 
fect on the pasta phase structure and on the size of the 
clusters. The largest effect occurs for the lower densities, 
namely the low density limit of the onset of the pasta 
phase, when a gas of free a-particles is considered: this 
defines an upper limit of the possible effect. The density 
at the upper border may be slightly larger for temper- 
atures and proton fractions close to the critical values 
above which the pasta phase disappears, like the results 
obtained for T=8 MeV. The density at the lower border 
is generally larger due to the softer EOS for the homoge- 
neous matter with a-particles. Including the interaction 
of the a-particles with the ui meson defines the exten- 
sion of the pasta phase to the limits obtained without 
a-particles. The effect of the model is mainly noticeable 
close the critical values of temperature and proton frac- 
tion: in particular for matter with trapped neutrinos at 
T=8 MeV we predict a pasta phase within TW, with or 
without a-particles, however, for NL3 there is a pasta 
phase only in the presence of non interacting a particles. 
A more precise knowledge of the effect of a-particles re- 
quires a self-consistent calculation, within, for instance a 
Thomas Fermi calculation. 

We have included in Tables HT1 and HTT1 the pressure at 
the transition between the pasta phase and the homoge- 
neous matter because it as been shown in [43| that this 
pressure defines the mass and moment of inertia of the 
crust of compact stars. The presence of a particles does 
not seem to have a large effect on the pressure at the 
transition. 

At this point a comment on the effect of thermal fluc- 
tuations is in order. The problem of the effect of ther- 
mal fluctuations on the pasta structures has been stud- 
ied by [H, |45| and it was shown that thermally induced 
displacements of the rod-like and slab-like nuclei could 
melt the lattice structure when these displacement are 
larger than the space available between the cluster and 
the boundary of the Wigner-Seitz cell. Using the elastic 
constants calculated in [44] the authors of [45] have calcu- 
lated the critical temperatures above which the ordered 
configuration is distroyed. While for the rod like clusters 
and for the densities and temperatures considered in our 
work the lattice would not be affected by the thermal 
fluctuations, for the slabs and densities considered T=5 
MeV is a limiting temperature. 

Table HVl allows a comparison between the models NL3, 
GM1 and TW at the crust-core transition for T=5 MeV 
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TABLE III: Densities of the inner edge of the crust (crossing points) for TW 



EOS 



p\ I pi P 1 I P 2 



pasta structure 



(fn 



(MeV / fm 3 



T = 5 MeV 
Y p = 0.5 



Y v = 0.3 



Y, = 0.4 



no as (*) 0.004 / 0.079 0.037 / 2.060 droplets, rods, slabs 

T„ a = 0.012 / 0.079 0.160 / 2.067 droplets, rods, slabs 

T„ a = 4r„ 0.004 / 0.079 0.038 / 2.060 droplets, rods, slabs 

no as (*) 0.003 / 0.086 0.018 / 1.209 droplets, rods, slabs 

I\, a = 0.007 / 0.086 0.049 / 1.211 droplets, rods, slabs 

T„ a = 4r„ 0.003 / 0.086 0.018 / 1.209 droplets, rods,slabs 

no as (*) 0.004 / 0.085 0.028 / 1.536 droplets, rods, slabs 

T va = 0.008 / 0.085 0.072 / 1.534 droplets, rods, slabs 

T va = 4r„ 0.004 / 0.085 0.026 / 1.537 droplets, rods,slabs 



T = 


= 8 MeV 












Y P 


= 0.5 


no as (*) 


"/- 


-/- 










r„„ = o 


0.030 / 0.042 


0.605 / 0.919 


droplets, 
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r — 4r 
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-/- 


-/- 






Y P 


= 0.3 
















no as (*) 


0.021 / 0.054 


0.221 / 0.714 


droplets, 


rods 






T va =0 


0.024 / 0.056 


0.263 / 0.747 


droplets, 


rods 






± va — 41 v 


0.021 / 0.054 


0.222 / 0.715 


droplets, 


rods 


Yf 


= 0.4 
















no as (*) 


0.023 / 0.050 


0.306 / 0.816 


droplets, 


rods 






Tva =0 


0.026 / 0.050 


0.389 / 0.818 


droplets, 


rods 






r — av 

J- va — ^t- 1 - v 


0.023 / 0.050 


0.306 / 0.816 


droplets, 


rods 



and several proton fractions. The values in this table in- 
clude the effect of interacting cv-particles. GM1 predicts 
the larger transition densities at the crust-core bound- 
ary. This reflects the larger binodal region this model 
presents. 



V. CONCLUSIONS 

In the present work we have studied the effect of a 
particles on warm low density stellar matter as found in 
the inner edge of the crust of a protoneutron star. We 
have considered three different types of relativistic nu- 
clear models: the parametrization NL3 [Tl| and GM1 [13[ 
of the non-linear Walecka model with constant couplings 
and the TW parametrization of the density dependent 
relativistic hadronic model [14| with density dependent 
coupling parameters. 

We have first considered an homogeneous neutral gas 
formed by protons, neutrons, electrons and a particles 
with a fixed proton fraction or in /3-equilibrium with 



trapped neutrinos. The a particles were described as 
a gas of bosons and two opposite situations were consid- 
ered: a free gas of a particles with binding energy 28.3 
MeV was used with no interaction taken into account, 
and a gas of particles with binding energy 28.3 MeV in- 
teracting through the exchange of uj meson which reduces 
the interaction to a repulsive interaction. It is expected 
that a realistic situation lies between these two extrema. 
The binding energy is density and temperature depen- 
dent and in the approach made in [241 ] the inclusion of 
a term that describes the temperature and density de- 
pendence of the binding energy is essential to dissolve 
the clusters. In the present work this was not consid- 
ered explicitly, because, in order to reduce the unknown 
a particle couplings, the attractive a meson-a interac- 
tion was not included. It was shown that the inclusion 
of the repulsive interaction is essential to avoid an over- 
prediction of a particles above p ~ 0.001 fm~ 3 and it 
is also the mechanism responsible for the dissolution of 
the a particle clusters in the present approach. In the 
low density limit (p < O.OOOlfm -3 ) all models considered 
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behave in a similar way: this behavior essentially coin- 
cides with the one predicted by the virial expansion [23[ 
when the interparticle interaction is negligible, and the 
system behavior is model independent. The main differ- 
ences between models occur for p > 0.001 fm~ 3 when the 
a particle fractions differences between models may be as 
large as one order of magnitude or even larger. TW pre- 
dicts the smallest fractions while GM1 the largest ones. 
Temperature shifts the maximum on the a particle dis- 
tribution and the density of cluster dissolution to larger 
densities, although the maximum values of the distribu- 
tions themselves decrease with temperature. The maxi- 
mum values of the a particle distributions also decrease 
when the proton fraction decreases. However, the proton 
fraction has no effect on the density localization of the 
maximum and neither on the density of dissolution of the 
clusters. 

We have next investigated the pasta phase with a par- 
ticles. We have described the pasta phase using the 
method described in [a, H B : the coexisting phases are 
determined from Gibbs conditions and surface energy 
and Coulomb interaction are added a posteriori. For the 
surface energy we have used the parametrizations deter- 
mined in terms of the proton fraction and the tempera- 
ture using a Thomas-Fermi calculation. It was shown in 
@, H @| that the method of the coexisting phases is very 
sensitive to the surface energy. We have performed our 
calculations for two temperatures T = 5 and 8 MeV, two 
proton fraction Y p = 0.5 and 0.3 and for /3-equilibrium 
matter with trapped neutrinos and a lepton fraction of 
0.4. We have analyzed the fraction of a particles as a 
function of density in the pasta phase and, comparing 
with the fraction of a particles in homogeneous matter, 
it was shown that it was larger by many orders of mag- 
nitude for densities above 0.01 fm~ 3 . It is important to 
stress that the prediction obtained within a homogeneous 
EOS calculation underestimates the a particle fraction. 
This certainly affects transport properties such as heat 
conduction and electrical conductivity. 

It was also seen that while for symmetric matter the 
a particle fraction decreases with temperature when in- 
teracting a particles are considered, for Y p — 0.3 the 



opposite occurs. This is an interesting effect related to 
the proton fraction in the background gas which increases 
with temperature for asymmetric matter. 

Finally we have analyzed the effect of the a particles on 
the pasta structure. It was shown that the effect is small 
except close to the critical temperatures and / or proton 
fractions when it may still predict a pasta phase while no 
pasta phase would occur in the absence of light clusters. 
A self-consistent calculation is necessary to give more 
quantitative predictions. Other small clusters should also 
be included in the calculation. 

VI. APPENDIX 

The surface tension coefficient, er, was calculated in 
the Thomas-Fermi approximation and fitted to the func- 
tional form given in, eq. (|41[) . with a{x)—<r{x,T = 0), 
a(T) and b(T) having the following expressions: 

Cr(x) — (To ex P( — °~\ £ 5 )(1 + 0,1 X + 0,2 x + 
0,3 X + 04 X + 05 X + Oq X ) , 



a{T) = a + ai T + a 2 T 2 + a 3 T 3 , 



b(T) = a Q + ai T + a 2 T 2 + a 3 T 3 . 

The results of the model fitting using the Thomas- Fermi 
approximation for the GM1, NL3 and TW parametriza- 
tions are given in table [V] for temperatures up to T = 10 
MeV. 
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TABLE IV: Densities of the inner edge of the crust (crossing points) for the T=5 MeV and r Q = 4T V 



Pi I A 

(fm- 3 ) 



pi i p2 
(MeV fm" 3 ) 



NL3 



GM1 



TW 



NL3 



GM1 



TW 



Y v = 0.5 
Y p = 0.3 
Yi = 0.4 



0.004 / 0.077 0.005 / 0.080 0.004 / 0.079 
0.003 / 0.082 0.004 / 0.088 0.003 / 0.086 
0.003 / 0.080 0.003 / 0.094 0.004 / 0.085 



0.039 / 2.065 0.049 / 2.076 0.038 / 2.060 
0.017 / 1.136 0.021 / 1.243 0.018 / 1.209 
0.024 / 1.430 0.028 / 1.527 0.026 / 1.537 



TABLE V: Surface tension coefficient parameters fitted within 
the Thomas-Fermi approximation for GMipjJ], NL3[11|] and 
TW[14|] parametrizations. 



GM1 a(x) 



a(T) 



b(T) 



(TO 
&l 

0<> 
Ol 
"2 
«3 

O4 

a 5 

(hi 



1.48801 
10.6647 

-4.32419 
61.6172 
-488.383 
1678.95 
-2450.41 
1296.86 



0.00521594 0.00640879 

0.0125187 0.000140886 

-0.00156999 -4.17615e-5 

8.4792e-5 1.6856e-6 



NL3 a(x) 



a(T) 



b(T) 



(TO 

(J\ 
Oo 

Ol 
2 

az 

4 

a 5 

Ofi 



1.12307 
20.7779 

-5.84915 
138.839 
-1631.42 
8900.34 
-21592.3 
20858.6 



0.0121222 0.00792168 

0.01664 -8.2504e-5 

-0.00137266 -4.59336e-6 

4.0257e-5 -2.81679e-7 



TW a(x) 



a(T) 



b(T) 



(Jo 
(T\ 

0() 
Ol 
2 

az 

4 

a 5 

0,6 



1.21736 
8.33425 

-2.18075 
18.0584 
-96.548 
259.517 
-296.69 
125.94 



-0.00481823 0.00844664 

0.0066498 -7.23379e-4 

-0.000267288 9.6817e-5 

1.34544e-5 -5.08488e-6 



